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In this note, we will discuss how to compactly express and prove the Jarzynski identity for an 
open quantum system with dissipative dynamics. We will avoid explicitly measuring the work 
directly, which is tantamount to continuously monitoring the system, and instead measure the 
heat flow from the environment. We represent the measurement of heat flow with Hermitian map 
superoperators that act on the system density matrix. Hermitian maps provide a convenient and 
compact representation of sequential measurement and correlation functions. 

PACS numbers: 05.30.Ch, 05.70.Ln 



Introduction 

When a classical system in thermal equilibrium is 
driven from that equilibrium by an external perturba- 
tion, then the irreversible work of that process is related 
to the system's free energy change by Jarzynski's equal- 
ity 



(e-^*^) = / p{W) e-P'^dW = e'P'^^. 



(1) 



Here, p{W) is the probability distribution of work W 
done on the system, (3 = I/ZcbT is the inverse tem- 
perature T of the environment in natural units, (fce is 
Boltzmann's constant) and AF is the Helmholtz free en- 
ergy change of the system. In other words, a Boltzmann 
weighted average of the irreversible work recovers the 
equilibrium free energy difference from a non-equilibrium 
transformation. 

The generalization of the Jarzynski identity to quan- 
tum dynamics |, i S i Et E III, B III [ll 113 , 

although technically straightforward, is not without some 
interesting conceptual complications. The total change in 
energy of the system /S.E = Q -\-W is equal to the work 
W applied via the time dependent perturbation plus the 
flow of heat Q from the environment. In a classical sys- 
tem, we can continuously measure the energy of the sys- 
tem and the ebb and flow of heat and work. However, 
continuous measurement is not possible in quantum me- 
chanics without severally disturbing the dynamics of the 
system. We should not attempt to continuously mon- 
itor the system [31, directly measure the work [I^ or 
radically redeflne work and heat [l^, [H, [23| • 

The simplest solution to this apparent dilemma is to 
disallow heat flow during the experiment [1, @ . The sys- 
tem is initially in thermal equilibrium with the environ- 
ment, but is decoupled and isolated from the environ- 
ment during the perturbation of the system. The work is 
then the difference in energy of the system between the 
beginning and end of the experiment. 



A more satisfactory solution is to realize that although 
we cannot continuously measure the work or system en- 
ergy, we can measure the heat flow from the environ- 
ment 0. The total system is divided into a system-of- 
interest and a thermal bath. This environment is large, 
rapidly decoheres and always remains at thermal equi- 
librium, uncorrelated and unentangled with the system. 
Consequently, we can measure the change in energy of the 
bath (i.e. —Q) without further disturbing the dynamics 
of the system. Essentially, we reexpress the open-system 
quantum Jarzynski identity as 



-/3AF 



(2) 



In this note I will discuss how to represent the heat flow 
and quantum Jarzynski identity with Hermitian maps, 
generalized measurement superoperators that act on the 
system density matrix alone. 



Background: Quantum Dynamics of Open Systems 

The mixed state of a quantum system that interacts 
with the environment can be described by a density ma- 
trix, p, a positive operator with unit trace. For an equi- 
librated system with Hamiltonian H interacting weakly 
with a thermal bath of temperature T the equilibrium 
density matrix is 



Pcq = 



(3) 



ilnZ 



where (3 — l/k^T, is Boltzman's constant 
trexp(— /3i7) is the partition function, and F = - 
is the Helmholtz free energy of the system. 

The dynamics of an open quantum system can be 
described by a quantum operation p' — Sp, a lin- 
ear, trace preserving, complete positive map of opera- 
tors [ll], [22, m, [13] . Any such complete, positive super- 
operator has operator-sum (or Kraus) representations. 



(4) 
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Conversely, any operator-sum represents a complete, pos- 
itive superoperator. The collection {^a} are known as 
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Kraus operators. The superoperator is trace preserving, 
and therefore conserves probabihty, if X]q ^q^" — I, 
where / is the identity operator. In the simplest case, 
the dynamics of a isolated quantum system is described 
by a single unitary operator = U^^. 

We are interested in the dynamics of a quantum system 
with a time dependent Hamiltonian, weakly coupled to 
an extended, thermal environment. Let the total Hamil- 
tonian of the combined system be 



eH" 



(5) 



where and are system and bath identity operators, 

is the time dependent Hamiltonian of the system, 
is the bath Hamiltonian, i?'"* is the bath-system inter- 
action Hamiltonian and e is a small coupling constant. 

We assume that initially the system and environment 
are uncorrelated and that the initial combined state is 
p|?q, where p^^ is the thermal equilibrium density 
matrix of the bath. We can derive a quantum operator 
description of the system dynamics by following the uni- 
tary dynamics of the combined total system for a finite 
time and then measuring the final state of the environ- 
ment. 



-T.ih\u[p' 




-l3Ef 



-{bf\U\h) p^ ihlU^bf) 



Here, U is the unitary evolution operator of the total 
system 



U : 



(7) 



tre is a partial trace over the bath degrees of freedom, 
{Ef} are the energy eigenvalues and {\bi)} are the or- 
thonormal energy eigenvectors of the bath, and is 
the bath partition function. For simplicity, and without 
loss of generality, we assume that the bath energy states 
are non-degenerate. 

It follows from the last line of Eq. ^ that the Kraus 
operators for this dynamics are 



A, 



-{b,\up 



(8) 



We further assume that the environment is large, with 
a characteristic relaxation time short compared with the 
relevant bath-system interactions and that the system- 
bath coupling e is small. Thus the environment remains 
very near thermal equilibrium, unentangled and uncorre- 
lated with the system. Consequently, the system dynam- 
ics of each consecutive time interval can be described by 



superoperators derived as in Eq. ([6]), which can then be 
chained together to generate a quantum Markov chain. 

p{t) = 5(t - 1, i) • • • S{s -I- 1, s -I- 2) S{s, s + 1) p{s) (9) 

In the limit of small time interval we obtain a continuous 
time quantum Markovian dynamics. 



p{t)^expi^J C{r)dTjpis) (10) 
where C is the Lindbladian superoperator [241] . 

Hermitian Maps and Sequential Measurements 

A measurement of a quantum system can be charac- 
terized by a collection of measurement operators {A^}, 
where A^Aa — I, and associated real valued mea- 
surement results, Ua- For example, the Hermitian op- 
erator H = of a standard von Neumann measure- 
ment can be decomposed into a collection of eigenvalues 
h and orthonormal projection operators P/i, such that 
H — ^i^hPfi- More generally, the measurement oper- 
ators of a POVM (Positive Operator Valued Measure) 
need not be projectors nor orthonormal [23 |. 

The probability of observing the ath outcome is 



Pa 



tr A^ypAl 



(11) 



and the state of the system after this particular interac- 
tion is 



/ 

Pa 



AapAl 



(12) 



The overall effect of the dynamics, averaging over differ- 
ent interactions, is the full quantum operation, Eq. (jl]). 

Rather than simple representing the effect of the mea- 
surement with the appropriate quantum operation, we 
can represent the effect and result of the measurement 
using a Hermitian map superoperator A : 



Ap = ^ UaAapAl^ 



(13) 



Note that this operator- value-sum cannot, in general, be 
recast as an operator-sum, since the measurement val- 
ues {tta} may be negative. An operator-value-sum maps 
Hermitian operators to Hermitian operators {H = H^) , 



[AH]^ - [J2 aaAa^HAl]^ - AH^ = AH 



(14) 



Conversely, any Hermitian map has an operator-value- 
sum representation 25]. 

Hermitian maps provide a particularly compact and 
convenient representation of sequential measurements 
and correlation functions. Let the Hermitian map A rep- 
resenting a measurement at time 0, B a different measure- 
ment of the same system at time i, and the quantum op- 
eration St represent the system evolution between these 
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two measurements. The expectation value of a single 
measurement is 

(a) ^tiAp^y^^ag tr AgpA]^ = p{a)aa (15) 



and the correlation function (fe(i)a(O)) can be expressed 
as 

{b{t)a{0)) = trBStAp{0) 

^Y^agbptvEfi [St[Agp{Q)Ai]] bI 

af3 

^'^p{a,P)aabi3 ■ (16) 

a/3 

(Note that expressions such as trABp, where A and B 
are Hermitian operators, often appear in perturbation ex- 
pansions and are frequently referred to as quantum cor- 
relation functions. However, since AB is not in general 
Hermitian these expressions do not directly represent a 
physical measurement.) 

Here we establish that, just as every Hermitian oper- 
ator represents some measurement on the Hilbert space 
of pure states, every Hermitian map can be associated 
with some measurement on the Liouville space of mixed 
states. Suppose that we have already decomposed the 
Hermitian map A into a value-operator-sum with val- 
ues {tta} and operators {Aq,}. Probability conservation 
requires that J2a ^a^a — I- If this condition is not 
met we can supplement the Kraus operators with an 
additional operator whose corresponding measurement 
value is zero. Note that Al^Aa is a positive operator 
and consequently B'^B = I — - ^a^a is also a pos- 
itive operator provided that c is a real number larger 
than the largest eigenvalue of J2a^a^a- Therefore, we 
can rescale the measurement outcomes {cAq} and Kraus 
operators {^}, append the additional operator B with 
measurement outcome 0, and associate the superoperator 
A with the measurement 



Ap^O BpB"! + 



Ac. Ai 

CUa —j= p — p 



(17) 



Note that the decomposition of a Hermitian map into an 
operator- value-sum representation is not unique [U [25| . 



Measurements of Heat Flow 

We can now construct a Hermitian map representa- 
tion of heat flow, under the assumptions that the bath 
and system Hamiltonians are constant during the mea- 
surement procedure and that the bath-system coupling is 
small. We construct a measurement on the total system 
and then project out the bath degrees of freedom, leaving 
a Hermitian map superoperator that acts on the system 
density matrix alone. 



The full, explicit measurement is 

h g f o d 

c b a be 

A 



We start with a composite system consisting of the bath, 
initially in thermal equilibrium, weakly coupled to the 
system (a). We measure the initial energy eigenstate of 
the bath (b), allow the system and bath to evolve to- 
gether for some time (c), and then measure the final 
energy eigenstate of the bath (d). The trace over the 
bath degrees of freedom (e) yields the final, unnormal- 
ized system density matrix, whose trace in turn (f) gives 
the probability of observing the given initial and final 
bath energy eigenstates. We then multiple by the Boltz- 
mann weighted heat (g) and sum over all initial and final 
bath states (h) to obtain the desired average Boltzmann 
weighted heat flow. 

The sum over initial states can be split into separate 
sums on the left and right projectors, since the bath is 
initially diagonal. Similarly, the sum over the final states 
can be split into separate sums on the right and left due 
to the final trace over bath degrees of freedom. As a con- 
sequence, we can rewrite the previous expression using 
the bath Hamiltonian. 



= trs trB[/^ 



We replace the bath Hamiltonian by = H — 
H^{t) (g) - ei/'"* [Eq. ©]. The total Hamiltonian 
commutes with the unitary dynamics and cancels, and 
we can discard the interaction Hamiltonian in the small 
coupling limit. 

« trs trB[e+f ® l""] U [e-^"" ® l""] [p^ ® • 



e 2 



e ' 2 



We now collect terms acting on the bath or system alone 



trs 



tre U [[e 



p^ e 2« ] 



and recover the Kraus operators {^q} describing the re- 
duce dynamics of the system, Eq. ([5]). 



trs J2 



e 2 



Finally, we have found that the average Boltzmann 
weighted heat flow can be represented by 



(18) 
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where S is the reduced dynamics of the system and the 
Hermitian map measurement superoperator TZ is 



Tltp 



P e 



(19) 



The paired Hermitian map superoperators act at the 
beginning and end of a time interval and measure the 
change in system energy during that intervaL This does 
not disturb the system beyond the disturbance aheady 
induced by couphng the system to the environment. 



Quantum Jarzynski Identity 

We are now in a position to derive the quantum Jarzyn- 
ski identity [Eq. ^] using Hermitian maps and the quan- 
tum operator formahsm. We spHt the total experimental 
time into a series of discrete intervals, labeled by the in- 
teger t. The system Hamiltonian is fixed within each 
time interval, and only changes in discrete jumps at the 
interval boundaries We can therefore measure the 
heat flow by wrapping the superoperator time evolution 
of each time interval St with the corresponding Hermitian 
map measurements TZ^^StTZt- Similarly, we can repre- 
sent the measurement of the Boltzmann weighted change 
in energy of the system with (exp(— /3A£') = trTZrSTZ^^ . 

The average Boltzmann weighted work of a driven, dis- 
sipative quantum system can therefore be compactly ex- 
pressed as 

(e-^^^) = tr [Ur (^[nt'StUt]^ 7^„- Vo'] (20) 
where p^'^ is the system equilibrium density matrix with 



system Hamiltonian Hf. 

Due to the structure of the energy change Hermitian 
maps TZ [Eq. (jl9p ] and the equilibrium density matrix 
[Eq. ^] this product of superoperators telescopes, leav- 
ing only the free energy difference between the initial and 
final equilibrium ensembles. 

(e-^'^) (21) 

= tr[7^,[7^;l5,7^,] • • • [7^^l527^2][7^^'5l7^l]7^^VS''] 

= tliTZriTZ-^SrUr] ■ ■ ■ [7^2"l527^2] [7^^'5l7^l] ^] 



trlTZrin-^Srllr] ■ ■ ■ [7^^ l527^2]7^^ '^i p^' 



z{oy 



= tr[7^,[7^-l5,7^,] • • • [7^-l527^2]7^^V^f|^] 



Zjr) 
Z{0) 



= exp{-/3AF} 



In the limit of small time intervals, we can express 
the quantum Jarzynski identity in the continuous time 
Llindblad form. 



(e-^^) = tr7^(^) exp <^ / R{t)-^ C{t)'R{t)dt \ 7^(0)-Vo'' 



= e-'^^^ (22) 
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